%% 10kV星接级联H桥SVG储能系统 - 初始化脚本
% 系统参数初始化
clear all; clc;

%% 1. 系统基本参数
% 电网参数
Vgrid_LL = 10000;                    % 电网线电压有效值 (V)
Vgrid_LN = Vgrid_LL/sqrt(3);         % 电网相电压有效值 (V)
f_grid = 50;                         % 电网频率 (Hz)
w_grid = 2*pi*f_grid;                % 角频率 (rad/s)

% SVG额定参数
S_rated = 12e6;                      % 额定视在功率 (VA)
Q_rated = 12e6;                      % 额定无功功率 (Var)
P_rated = 5e6;                       % 额定有功功率 (W)
P_storage_ratio = 0.4;               % 储能功率占比

%% 2. 主电路参数
% H桥级联参数
N_bridge = 11;                       % 每相H桥个数
Vdc_ref = 950;                       % 单个H桥直流侧电压 (V)
Vdc_total = N_bridge * Vdc_ref;     % 每相总直流电压 (V)
N_levels = 2*N_bridge + 1;           % 电平数 (23电平)

% IGBT参数
V_IGBT = 1700;                       % IGBT额定电压 (V)
fsw = 2000;                          % 开关频率 (Hz)
Tsw = 1/fsw;                         % 开关周期 (s)

% 滤波参数
L_filter = 0.15;                     % 滤波电感 (mH/MVA基准)
L_ac = L_filter * (Vgrid_LN^2/S_rated) * 1e-3;  % 实际滤波电感 (H)
R_ac = 0.01 * (Vgrid_LN^2/S_rated);  % 滤波电阻 (Ω)

%% 3. 储能系统参数
% 超级电容参数
Vsc_rated = 1920;                    % 超级电容额定电压 (V)
C_sc = 15;                           % 超级电容容量 (F)
SOC_init = 0.7;                      % 初始SOC
SOC_max = 0.9;                       % 最大SOC
SOC_min = 0.2;                       % 最小SOC

% DC/DC变换器参数
eta_dcdc = 0.95;                     % DC/DC效率
P_dcdc_max = P_rated/3/N_bridge;     % 单个DC/DC最大功率 (W)

%% 4. 控制系统参数
% 电流环PI参数
Kp_i = 0.5;
Ki_i = 50;

% 电压环PI参数
Kp_v = 2;
Ki_v = 20;

% 直流电压控制PI参数
Kp_dc = 0.1;
Ki_dc = 10;

% SOC均衡控制参数
Kp_soc = 0.05;
Ki_soc = 0.5;

% PLL参数
Kp_pll = 180;
Ki_pll = 3200;

%% 5. 负载参数
% 总负载
P_load_total = 3.5e6;                % 总有功负载 (W)
Q_load_ind = 1.25e6;                 % 感性无功 (Var)
Q_load_cap = 0.25e6;                 % 容性无功 (Var)
Q_load_total = Q_load_ind - Q_load_cap;  % 净无功负载 (Var)

% 负载分配
P_industrial = 0.6 * P_load_total;   % 工业负载 60%
P_commercial = 0.26 * P_load_total;  % 商业负载 26%
P_residential = 0.14 * P_load_total; % 居民负载 14%

% 负载阻抗计算
Z_base = Vgrid_LN^2 / (S_rated/3);   % 基准阻抗
R_load = Vgrid_LN^2 / P_load_total * 3;  % 等效电阻
X_load = Vgrid_LN^2 / Q_load_total * 3;  % 等效电抗
L_load = X_load / w_grid;            % 等效电感

%% 6. 电网强度参数
% 强电网
SCR_strong = 10;                     % 强电网短路比
X_grid_strong = 1/SCR_strong;        % 强电网阻抗 (pu)
L_grid_strong = X_grid_strong * Z_base / w_grid;  % 强电网电感 (H)

% 弱电网
SCR_weak = 2;                        % 弱电网短路比
X_grid_weak = 1/SCR_weak;            % 弱电网阻抗 (pu)
L_grid_weak = X_grid_weak * Z_base / w_grid;      % 弱电网电感 (H)

% SCR切换阈值
SCR_threshold = 3;                   % 强弱电网切换阈值

%% 7. 仿真参数
Ts = 1e-5;                          % 仿真步长 (s)
T_sim = 3;                          % 仿真时间 (s)

% 事件时间设置
t_switch = 0.5;                     % 跟网/构网切换时间 (s)
t_disturbance = 1.0;                % 扰动发生时间 (s)
t_recovery = 2.0;                   % 扰动恢复时间 (s)

%% 8. 测量和显示参数
% THD限值
THD_limit = 5;                      % THD限值 (%)

% 频率误差阈值
f_threshold = 0.1;                  % 频率误差阈值 (Hz)

% 功率基准值
S_base = S_rated;                   % 视在功率基准
V_base = Vgrid_LN;                  % 电压基准
I_base = S_base/(3*V_base);         % 电流基准

%% 9. 初始化控制变量
% 控制模式 (0=跟网, 1=构网)
control_mode_init = 0;

% 参考值初始化
Vdc_ref_init = Vdc_ref * ones(3, N_bridge);  % 初始直流电压参考
P_ref_init = -P_rated * 0.2;        % 初始有功参考 (充电)
Q_ref_init = Q_load_total * 0.8;    % 初始无功参考

%% 10. 保存参数
save('SVG_Parameters.mat');
disp('系统参数初始化完成！');